UNIVERSITATEA TEHNICA DE CONSTRUCTII BUCURESTI 


Catedra de Hidraulica si Protectia Mediului 


Stratul limita laminar si turbulent 
dezvoltat la curgerile exterioare 


CONDUCATOR DOCTORAT 
Prof. Dr. Ing. Mircea Degeratu 


DOCTORAND 
Ing. Tudor Baracu 


2011 


Cuprins 


Introducere si istoric al teoriei stratului limita 

Ecuatiile generale ale stratului limita stationar 

Stratul limita pe o suprafata 

Ecuatiile stratului limita de-a lungul unei placi plane 

Grosimea stratului limita 

Investigatii experimentale ale stratului limita 

Forma adimensionala a ecuatiilor Navier-Stokes la curgerea plana 
Curgerea potentiala si stratul limita in apropierea bordului de atac al placii 
Stratul limita turbulent 

Stratul limita laminar turbulent cu gradient de presiune 

Ecuatiile integralei impulsului si ale integralei energiei pentru stratul limita 
Desprinderea stratului limita la curgerea in jurul unui corp 

Exemplu de simulare numerica a curgerii bidimensionale in jurul unei elipse 
Bibliografie 


Introducere si istoric al teoriei stratului limita 


Prandtl in 1904 a introdus teoria stratului limita cu scopul de a depasi dificultatea rezolvarii 
ecuatiilor ccurgerii potentiale, cand numarul Reynolds este foarte mare. Aceste ecuatii rezultate prin 
neglijarea vascozitatii nu satisfac conditiile la limita pe o suprafata solida. Reconsiderand curentii 
irotationali, Prandtl a postulat existenta intre ei si suprafata solda a unui strat subtire de fluid, in care 
viteza creste foarte repede de la zero la perete, la valoarea teoretica din curent. Intr-un astfel de strat in 
care derivata spatiala a vitezei este mare, fortele de vascozitate sunt foarte mari. Lucrarile 
experimentale au confirmat corectitudinea ipotezelor lui Prandtl aratand nu numai ca stratul limita 
exista, dar si ca distributia vitezelor prin el se dezvolta in concordanta cu mecanica dezvoltata ulterior. 
De asemenea, aceste experiente au mai aratat si in ce conditii stratul limita dezvolta instabilitatii si 
devine turbulent. Prandtl (1914, 1927) s-a dovedit a fi si descoperitorul efectelor acestor tranzitii 
asupra mecanicii stratului (impiedicarea desprinderii) si asupra variatiei cu numarul Reynolds a 
diferitelor modele de curent din jurul corpurilor. 

Ipoteza stratului limita ne ajuta sa reconciliem intuitia ca efectele vascozitatii asupra stratului limita 
sunt neimportante cand vascozitatea este mica, cu conditia de nealunecare care trebuie sa fie 
satisfacuta pe frontiera solida, oricat de mica ar fi vascozitatea. Aceasta reconciliere a fost obiectivul 
principal al lui Prandtl si a constituit piatra de hotar in dezvoltarea mecanicii fluidelor. Stratul limita 
este in definitiv stratul in care viteza fluidului face tranzitia la valoarea zero ceruta pe frontiera solida 
la o valoare finita pentru fluidul fara frecare. In apropierea frontierei corpului cresterea vitezei 
(tangentiale) este rapida, ceea ce conduce la o valoare ridicata a derivatei acesteia dupa normala la 
suprafata (gradientul). Din acest motiv si tensiunile tangentiale ating valori considerabile in stratul 
limita. Dimpotriva, in afara acestui strat derivata normala a vitezei are o valoare redusa si din aceasta 
cauza tensiunile tangentiale pot fi neglijate. Simplificarea introdusa in felul acesta in tratarea 
matematica a problemelor de miscare in care numarul Reynolds are valori mari este foarte importanta 
si numeroase probleme au putut fi rezolvate pe aceasta cale. 


Un scurt istoric al teoriei stratului limita: 

e In 1904 Prandtl prezinta lucrarea “Miscarea fluidelor cu vascozitate mica”. Aceasta lucrare a 
fost o scurta expunere a teoriei stratului limita, separarea acestuia, si intarzierea (controlul) 
separarii in stratul limita prin efectul indus de suctiune. 

e In 1908 Blasius prezinta rezultate ale stratului limita peste o placa plana intr-un curent uniform, 
si peste un cilindru. 

e In 1908 Bolze raporteaza rezultate ale stratului limita pentru o sfera. 

e In 1910 Prandtl aplica conceptul de strat limita problemelor de transfer de caldura 

e In 1914 Prandtl explica reducerea rezistentei prin frecare la curgere in jurul unei sfere dupa 
atingerea unei anumite valori pentru numarul Reynolds. 

e In 1914-1924 Von Karman dezvolta ecuatiile integrale ale impulsului de-a lungul stratului 
limita. 

e De-a lungul perioadei care a urmat ulterior pana in prezent au mai fost abordate: 

o Predictii ale tranzitiei sub influenta parametrilor de gradient al presiunii, curbura, 
compresibilitate, rugozitate 

o Stratul limita turbulent 

o Principii de masuratori ale stratului limita (prin anemometru cu fir incalzit, anemometru 
cu laser Doppler — LDA, velocimetre cu imagine de particule — PIV 

o Stratul limita tridimensional 


o Stratul limita nestationar 

o Stratul limita atmosferic 

o Stratul limita computational: Ecuatiile integrale ale impulsului, modele de turbulenta ca 
RANS, DNS, LES, reducerea frecarii la curgere, etc 


Ecuatiile generale ale stratului limita stationar 


Trebuie adagat ca ecuatiile stratului limita implica ecuatiile generale ce stau la baza curgerii 
fluidelor, la care se adauga conditia limita de viteza nula la perete. Comportamentul neliniar al miscarii 
fluidelor implica in stadiile ulterioare un aparat matematic complex pentru rezolvarea si descrierea el. 


Ecuatia continuitatii 


O(pu) 0(pv) lpw) _ i 
ax dy ðz 


Ecuatiile Navier-Stokes de miscare a fluidelor 
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Stratul limita pe o suprafata 


S-a constatat din o serie de solutii exacte ale ecuatiilor Navier-Stokes ca grosimea stratului limita 
este proportionala cu radacina patrata a vascozitatii cinematice: 
d~Vv 
Atata timp cat componenta vitezei paralela cu peretele creste de la valoarea zero pe perete la 
valoarea | in liniile de curent de-a lungul stratului de grosime 6, avem: 


Ou 1 

dy 6 
d*u 1 
d*y 62 


Daca aceste relationari sunt combinate cu ecuatiile Navier-Stokes se obtine in final estimarea 


grosimii stratului limita: 


Ecuatiile stratului limita de-a lungul unei placi plane 


Problema miscarii fluidului pe o suprafata plana a fost prima data abordata de H. Blausius. 
Reprezinta cel mai simplu caz de abordare pentru stratul limita, si in acelasi timp reprezinta un 
important caz practic. 


In cazul acestui curent uniform pe placa plana nu exista punct de desprindere deoarece dp/dx=0, 
astfel incat avem un strat limita pur. 


Ecuatia continuitatii 


Ecuatiile Navier -Stokes 


Conditiile limita 


In situatia unei accelerari bruste a unei placi s-a gasit empiric ca: 
6~VVT 
unde T este timpul in care s-a facut miscarea. 
Putem asimila pe T ca fiind timpul consumat pentru traversare de la marginea accelerata pana la 
distanta x de aceasta. Pentru o particula din afara stratului limita acesta este: 


=T 


astfel incat putem pune: 


Vx 
en ee 


unde Uœ reprezinta viteza curentului la mari distante. 


o 


Uz U= const? 


Frontiera” strotului limită 
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substratul lominor 
zona turbutenta 


zone laminară 
franzitie 


Dupa cum este reprezentat in imaginea de mai sus, pentru o anumita lungime curgerea este laminara 
prin dominarea fortelor vascoase asupra celor inertiale. De la un anumit punct incepe o zona de 
tranzitie cand stratul limita se separa. O ultima zona este de turbulenta, vectorul viteza incepe sa se 
uniformizeze de-a lungul sectiunii de curgere. 


Acum putem introduce coordonata adimensionala 
~Y 
a 
astfel incat 
Uss 
ȘI e PT 


Introducem functia de curent ¥ (x, y) deci: 


Y= J/vxUof(n) 


unde f (n) denota functia de curent adimensionala. Astfel, componentele devin: 


OY  0W On ; 
oy On ay Usf (n) 
_ Ov 1 na ; ) 
vi Ox 2 x MS 


Introducand u si v in ecuatia Navier-Stokes scrisa mai inainte obtinem: 
2 


U% Ign + Ves ( 1 ) n __ Ux mr 
snuff += (nf - Af" =v f 
Dupa inlocuiri si simplificari se obtine ecuatia lui Blausius: 


ff" + 2p" =0 


Grosimea stratului limita 


Tranzitia vitezei fluidului de la perete spre exterior se realizeaza in mod asimptotic si abrupt spre 
viteza straturilor de fluid exterioare de-a lungul unei distante mici de perete, deci in consecinta este o 
variatie neliniara a acesteia. Acest lucru determina o evaluare aproximativa si nu exacta a grosimii 
stratului limita. Astfel se considera o grosime a stratului limita arbitrara distanta fata de perete in care 
viteza fluidului atinge 99% din viteza exterioara. Uneori in locul notiunii de grosime de strat limita 
mai este utilizata notiunea “grosimea de deplasare” care o vom nota cu 6, . Aceasta este distanta de 
deplasare in care campul potential de curgere exterioara este deplasat in afara ca o consecinta a 
descresterii vitezei in stratul limita. 

Astfel grosimea de deplasare a stratului limita va fi caracterizata de ecuatia: 


6 = Ju —u)dy 
0 


Obtinerea exacta a grosimii stratului limita este un demers utopic intrucat sunt foarte multi parametri 
care o influenteaza, dintre care multi sunt haotici, iar in acelasi timp vascozitatea fluidului in stratul 
limitavariaza descrescator spre exterior si asimptotic. 

Astfel componenta paralela cu peretele u la o distanta mare de perete va tinde asimptotic catre 
viteza globala U» (in care functia f'(n) tinde asimptotic la 1. Se obtine o viteza u = 0.99 U,, pentru 


_ [Vă Sâ 
y= 4] aa TA 


Descresterea volumului de curgere datorita influentei frecarii este: 


n = 5 si atunci: 


|. u)dy 
0 
ca urmare vom avea relatia: 
iUs = Í (U» — u)dy 
0 


sau 


r u 
A co 
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Dar se cunoaste ca: 


z=) 


Us x = 
pe e [terete = fim- ran 
ome =0 


unde 7, denota un punct din exteriorul stratului limita.Din tabelul care exprima pe 


n, f(y), f'(m) si f'"(m)se poate obtine ca: 
6, = 1.7208 ead 
171. U, 


deci 
In figura de mai jos este prezentata distanta y = 6, in care liniile de curent externe ale potentialului de 


curgere sunt deplasate datorita efectului de frecare de langa perete. Grosimea stratului limita 6 va fi de 
trei ori mai mare decat grosimea de deplasare 6,, deci 6 = 3 ô. 
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Putem sa evaluam o a doua marime numita “grosimea de impuls” 6, . Pierderea de impuls in stratul 
limita comparat cu potentialul de curgere este data de valoarea: 


(oe) 


| u(U» — u)dy 


0 
Astfel relatia pentru aflarea impulsului corespunzator grosimii este: 


(0.0) 


p UZ, 62 =p i u(U,, — u)dy 
y=0 


CO 
pe | ae 
U, y 
y=0 
in urma substitutiilor se obtine: 


62 = P | rou — f'(n)] dn = maat 
n=0 
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Sau 


Trebuie mentionat ca teoria stratului limita se aplica doar pentru anumite intervale ale numarului 


Reynolds = == x 
Langa marginea frontala a placii teoria stratului limita nu se mai aplica intrucat nu mai are loc 
conditia: 

0*u 0*u 

9x2| ~ lay? 


si atunci se pot folosi calculatii elaborate de G. F. Carrier, C. C. Lin, B. A. Boley si M. B. Friedman. 


Investigatii experimentale ale stratului limita 


Un bogat material experimental elaborat in urma cercetarilorstratului limita a fost propus de J. 
Nikuradze. S-a demonstrat ca formarea stratului limita este foarte mult influentata de forma muchiei 
frontale a profilului si de asemenea de un foarte mic gradient de presiune care poate exista in curgerea 
externa. 

Pentru o serie de distante fata de varful profilului, Nikuradze a elaborat graficul de mai jos. 
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Se observa cum teoria lui Blausius este confirmata de rezultatele experimentale. 


Forma adimensionala a ecuatiilor Navier-Stokes la curgerea plana 


Teoria stratului limita poate fi vazuta in sensul integrarii asimptotice a ecuatiilor Navier-Stokes. 
Aceasta ipoteza ne conduce la ideea unei relationari intre numarul Reynolds si caracteristicile stratului 
limita din jurul unui corp in cauza. 

Dependenta caracteristicii stratului limita de numarul lui Reynolds se poate expune prin 
transformarea Ruark centro-afina. 

In ecuatiile Navier-Stokes si a continuitatii se vor folosi marimi adimensionale astfel: 


Obtinem la o curgere stationara bi-dimensionala: 


Ou’ j du dU, 1 8 u 
“Jx” dy’ — dx! Rey? 
du 9% 
ax’ ðy 
Cu conditiile limita: 
(e = 0 > u=v=0 
y' — © = u' = U'(x") 
UoolL 
In care Re = i 
Daca punem: 
Vv |UyL 
v" = vV Re = — |—— 
Uv v 
UL 
y" = y’ /Re — A et eee 
Ecuatiile Navier-Stokes si a continuitatii devin: 
ðu' „ðu dU’ 9u' 
© ax! dy" dx’ dy" 
du’ Ov" 
ax’ Oy 
Cu conditiile de granita: 
A 0 => uU' — y" — Q 
y" > u' = U'(x") 


Solutia la aceste ecuatii este valabila pentru orice numar Reynolds in conditiile de strat limita laminar. 


10 


Curgerea potentiala si stratul limita in apropierea bordului de atac al 
placii 


Discrepantele observate intre teorie si si rezultatele experimentale provin, pe de o parte, din faptul 
ca in calcule s-a presupus placa de lungime infinita, iar, pe de alta parte, din cauza ipotezei ca stratul 
limita incepe chiar de la bordul de atac al placii. In realitate, perturbatia produsa de placa se manifesta 
Si inaintea acestui bord, astfel incat studiul miscarii in vecinatatea lui x=0 cere un studiu special. Dintre 
cercetarile facute asupra miscarii in vecinatatea bordului de atac se amintesc acelea ale lui G. F. 
Carrier si Lin. 

Consideram ecuatiile stratului limita in curgerea plana stationara: 
ðu du dU du 


Ox dy dx dy? 
ðu i Ov i 
Ox Oy 


cu conditiile limita: 
f y=0 =a u=v=0 


y = 00 > u =U 
Ecuatia continuitatii este integrata prin introducerea functiei de curent ¥ (x, y) in care: 
OV OV 
u = ay ; PV = — Ae 
De asemenea: 
z=xt+iy=re 
oY 
w = u + iv = -2i — 
Oz 


Astfel ca ecuatia de miscare devine 
OV 92W away _ dU ð? y 


dy ðy? ax ðy? — aera dy? 


Cu conditiile limita 


OF OW 
y=0 > Di ay = 0 
OV 
y = 00 o ay = 
Vom face transformarile de variabile utilizand si un factor de scara g(x). 
x yvRe 
pe au) 
Functia de curent va deveni adimensionala prin substitutia: 
Gs y)VRe 
fn) = PU) a) 
Vom avea 
OV Of E 
u = ay =U an = Uf 
—VRev = VRe = Lf —(Ug) + Ug & Le nf’) 


Obtinem urmatoarea ecuatie diferentiala in coordonate adimensionale: 
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U Of! O 
pr apps pap = 2ar( rb pd) 


Unde a si P sunt contractiile pentru urmatoarele functii de x: 


_Lgd L dU 
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gdi U9)» p RIF 
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Cu conditiile de limita 
n= 0 > ff =0 
ate Ss fad 
Functia de curent adimensionala f (n) va trebui sa satisfaca ecuatia ordinara cu derivate partiale: 
f" + aff'+BU- f?) =0 
Cu conditiile de limita 
n=0 > f=f=0 
fon a fed 
Din expresiile lui a si 6 obtinem: 


L d 
2 a — B =——(g*U 
a—B U dx 9 ) 
si daca 2a — p + 0 atunci 
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Mai departe avem ca 
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ea Dati 2y 2 gf 
U U% g g 


Dupa integrare obtinem ca 


U 
e: = KaB 
(Z) g 
unde K este constant. 
Se obtine in final ca 
1 
xXx/UN\2 
- female) 
g (2a 6)? TA 
B 
U mee X12a—B 
ea 7 ap) — 
Us KF |Qa P7] 
Se face notatia 
£ 2m 
= m+1 
Atunci 
m 
KGtm) ( 2 =) 
= 1+mL 
2 xU, 
g= r 
1+mL U 


Astfel pentru coordonata n se obtine transformarea 
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_ m+1 U 
n=y 2 vx 


O solutie polinomiala satisfacand conditiile de mai sus pentru o placa plana este considerand 
Y = Axy? 
Vor apare regiunile ca in figura: 


- Regiunea I- a solutiei potentiale, unde y > 1 


- Regiunea II — in care y ~ 1 siy &x 

- Regiunea III — in care y « 1 

- Regiunea IV in care ¥ nu are o solutie analitica simpla 

Un caz particular este curgerea in jurul unei pene triunghiulare. Se ia: 
U(x) = Cx" 


2 


—— — M 


= 


Unde C este o constanta. 
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Stratul limita turbulent 


Tratarea stratului limita turbulent este de departe mai dificila datorita variatiilor dependente de timp a 
proprietatilor curgerii. Una din cele mai larg utilizate metode in abordarea curgerilor turbulente este 
aplicarea decompozitiei Reynolds in care marimile de curgere sunt descompuse in componenta medie 
si fluctuanta. 


Ecuatia continuitatii in stratul limita turbulent 


ða 9 ôw 


Ecuatiile Navier-Stokes in stratul limita turbulent (in forma RANSE) 


oa(u) au) oa(u) 1 0p _ 1fa =E 0 eS 0 -e 
— jaar se = E, ee A aan] (eee, cree Ta NE (SA Ta ees een Fart: 
Ax + Vv ay + WwW Az aT ama a+ [= ( puu Jez pv'u') += ( pw')| 


a(v) a) ð (T) 10p = ya coe 0 sn | SNI 
— a ke Ses — = E Na die. Se A a (sd ARE ral EA ral AS fa rar 
Ax v ay + WwW 37 ay" ek puv)+ al pvv)+ > pw'v)| 


0(w) o (w) o (w) 1 0p _ 19 == O == 0 == 
= pq a a — i pe aaa E A | ( A I — (eee A I ee pe 1 I 
u Ax +v 3y +w Az Da at [al puw)+ al pv'w') + = ( pw'w)| 


Termenii aditionali de tipul v'u' din ecuatiile stratului limita turbulent sunt cunoscuti si ca tensiuni 
taietoare Reynolds, si sunt necunoscute apriori. Solutia ecuatiilor stratului limita turbulent de aceea 
necessita folosirea unui model de turbulenta care are in vedere sa exprime tensiunea taietoare Reynolds 
in termeni de variabile de curgere cunoscute sau derivate. Lipsa de acuratete si generalitate ale acestor 
modele este un obstacol major in predictia cu succes a proprietatilor de curgere turbulenta in dinamica 
fluidelor moderna. 


Turbulenta este generata de fortele taietoare. Cu cat fortele taietoare sunt mai mari, cu atat 
turbulenta este mai puternica. In imediata apropiere a peretelui se formeaza un substrat laminar de 
grosime 6,. Cu cat nivelul general de turbulenta a curgerii este mai mare, cu atat grosimea acestui 
substrat scade. 


ð 


u/ U, 1 


Maximum de turbulenta se atinge in zona in care fortele taietoare sunt maxime. 
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Schema de la care se pleaca pentru studiul turbulentei este ca in figura de mai jos. 


Substratul laminar 


este in imediata vecinatate a peretelui 
regim laminar, straturile de curent sunt paralele 
efect predominant al vascozitatii asupra efectelor inertiale moleculare si de amestec 


Strat de tranzitie 


asigura trecerea de la substratul laminar la zona de turbulenta 

uneori acest strat este neglijat, considerandu-se o trecere directa de la substratul laminar la zona 
turbulenta, fapt preferat in modelarea matematica 

regimul este tranzitoriu, straturile de fluid incep sa se intersecteze sau indeparteze local 

schimb masic permanent si relativ stabil cu zona turbulenta 

efectele vascozitatii sunt aproximativ egale cu cele inertiale moleculare si de amestec 


Zona turbulenta 


straturile de fluid se intersecteaza in toata masa fluidului din aceasta zona, amestecul este 
accentuat, distributie vartejuri. 
Efectele vascozitatii sunt dominate de cele inertiale moleculare si de amestec 


Stratul limita laminar turbulent cu gradient de presiune 


Experimente mai timpurii in acest domeniu au fost facute pe pereti plani de catre Nikuradze, J. , 


Doench, F. , Hochschild, H. etc. Aceste experimente au demonstrat ca forma profilului de viteza 


depinde foarte mult de gradientul de presiune. 


Numere caracteristice 
Pentru a da o descriere a comportamentului esential al stratului de viteze este necesar sa cunoastem 


grosimea si sa avem un indiciu al distributiei vitezelor in stratul limita. 


Atata timp cat grosimea stratului limita d(x) de-a lungul caruia stratul disipativ se uneste cu 


curgerea externa fara frecari U(x), astfel incat u(x, y = 6) = U(x) nu se poate defini cu acuratete. De 


aceea se definesc 3 marimi diferite: 


6, =f N (1 — z) dy - grosimea de deplasare 


62 =] ° (1 — =) zdy - grosimea de impuls 


2 
= fe E — (=) |ay - grosimea de energie 
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Aceste cantitati pot fi facute adimensionale introducand numere Reynolds derivate din ele formate cu 
viteza externa. 


Re, =6,—; Rez = 63— 
2 27 3 a 
Profilul de viteze depinde puternic de gradientul extern de presiune, exprimat prin derivata dU/dx, 
si este caracterizat de un numar al factorilor de forma. Acestea sunt facute de asemenea 


adimensionalede preferinta si pot fi definite in forma de rapoarte ale grosimilor. Se obisnuieste sa se 
foloseasca contractiile: 


Masuratorile au indicat ca profilele turbulente ale vitezelor pot fi descrise aproximativ de catre o 
familie parametrica de curbe. 
E. Truckenbrodt a introdus factorul de forma modificat 


[H32 dH32 [H23 dH23 
(H12)00(H21—-1)H32 + e (H12)oo (H12—1)H32 


H =e 

Valoarea de referinta (H23).. = (1/H32).. a fost aleasa ca cea mai mica limita a integralei pentru 
ca reprezinta o valoare medie pentru curgerile fara gradient de presiune. In caz de strat limita turbulent 
se ia (Hi2)œ = 1.3. 
In caz de curgere cu gradient de presiune zero se gaseste ca H = H,, = 1 prin definitie (valoare medie 
in caz de curgere turbulenta). 
Curgerile cu gradient de presiune contrar (presiunea crescand in directia curentului aval) sunt 
caracterizate de H, < H < 1 in timp ce pentru curgerile accelerate (presiunea descrescand) se gaseste 
ca 1 < H < Ho, unde H, este factorul de forma pentru profilul de viteze cu separare incipienta, si Ho 
reprezinta factorul de forma a profilului de curgere stagnanta bi-dimensionala. 
Dupa K. Wieghardt, factorii de forma H42 si H32 sunt legati unul de altul prin ecuatia 


Hız = nE 
(3H32 — 4) 
cu asumarea numita legea puterii profilului. 
Valori numerice indicate in literatura pentru stratul limita turbulent pentru care separarea are loc 
variaza considerabil. J. C. Rotta recomanda 4.05 > (H,,), > 4.0 sau H, = 0.723. 
In figura de mai jos este prezentat raportul de grosime al stratului limita H32 = 63/62 versus H42 = 


6,/62 conform lui J. C, Rotta si W. Wieghardt 


a 


Rotta 
——— Wieghardt 
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Ecuatiile integralei impulsului si ale integralei energiei pentru stratul 
limita 


O calculare completa a stratului limita cu ajutorul ecuatiilor diferentiale este atat de complicata 
incat poate fi abordata doar prin intermediul tehnicilor computationale, dar si atunci poate fi mare 
consumatoare de timp pentru a obtine si o buna exactitate. 

Ecuatia ce va fi implicata va fi cea a integralei impulsului pentru stratul limita, sau ecuatia integrala a 
lui Karman. 


ee a . , r 
Efortul tangential la perete, Tọ poate fi substituit cu (u =) astfel incat va fi valida atat pentru 
0 


curgerea laminara cat si pentru cea turbulenta, cu conditia ca intr-un caz ulterior u si v vor fi raportate 
la timpul mediu al componentei de viteza respective. 


Componenta de viteza normala, v, poate fi inlocuita prin: 
y 


=~ | (az) 
a ax)“ 
0 
Astfel, vom avea: 
h y 
| ðu =f (Sa Pia Ie To 
“Jx dy) NOx y ax |“ 7 p 
y=0 0 
Integrand prin parti obtinem pentru termenul al doilea de sub integrala: 
h y h h 
f SE 4 =u f (S)e du 
Oy J NOx a a ax) 7 “Ix 
y=0 0 0 0 


astfel ca 


h 
fh ðu oo ua = To 

"Ox Ox Ox y= 
0 


care poate sa fie comprimata in: 


h h 

o dU To 
| putu -wlay += f U -wdy = D 
0 0 


Atata timp cand in ambele integrale integranzii se anuleaza in afara stratului limita, este permis sa 
punem h > œ. 
Introducem acum grosimile de deplasare 6, si grosimea de impuls 62. Ele sunt definite prin: 


6,U = fu —u)dy, 6,U = [uw —u)dy 
y=0 y=0 
In ecuatia integrala a lui Karman, integralele sunt indepente de x, cat si limita superioara h. Se obtine 
in final: 
T a dU 
— = — (U*6,) + ôU — 
reer (02) + ô1U e 


Aceasta este ecuatia impulsului-integral pentru stratul limita bi-dimensional, incompresibil. Atata timp 
cat nu este facut nici un statut anume pentru Tg, ecuatia este aplicabila atat stratului limita laminar cat 
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si turbulent. Aceasta forma a ecuatiei impulsului-integral a fost data de H. Gruschwitz. Ea isi gaseste 
aplicabilitatea in teoriile de aproximare a straturilor limita laminar si turbulent. 

Folosind o abordare similara, K. Weighardt a dedus o ecuatie integrala pentru stratul limita laminar. 
Aceasta ecuatie este obtinuta multiplicand ecuatia miscarii cu u si apoi integrand de la y=0 la 
y=h>5(x). Substituind din nou v din ecuatia continuitatii obtinem: 


y 
»| po us HSE geca jg a, 
Ax y u FE Yy=u|u y 
0 


Al doilea termen poate fi transformat integrand prin parti: 


h 
ðu 
De an 
fhal [5 dy ||dy=3 | u) a dy 
0 


Si in timp ce combinam termenul unu si trei obtinem 


h h 
dU 1 O 
Eanes | Aa Bad re ee e. 
fhs u’ = u dy = 5 | ut U“)dy 
0 0 
In final, integrand partea dreapta prin parti se obtine 
1 d r r ðu 
paan 
zP | WU —udy =p By 
0 0 


Limita superioara a integralei poate fi si aici de asemenea inlocuita cu y = ©, pentru ca integranzii 
devin egali cu zero in afara stratului limita. Cantitatea u(du/dy)* reprezinta energia pe unitatea de 


volum si timp, care este transformata in caldura prin frecare (disipare). Termenul 5 p(U 2 — u?) din 

partea stanga reprezinta pierderea in energie mecanica (cinetica si energie de presiune) avand loc in 
E 1 d 

stratul limita raportat la potentialul de curgere. In timp ce termenul i ies f a u(U* — u?)dy reprezinta 


fluxul de disipare a energiei, si partea stanga reprezinta rata de schimbare a fluxului energiei disipate 
pe unitatea de lungime in directia x. 

Daca in plus la grosimile de deplasare si de impuls introducem grosimea energiei de disipare, 63, prin 
definitie atunci avem: 


== i u(U* — u*)dy (grosimea de energie) 


Si atunci putem rescrie ecuatia integrala a energiei in forma simplificata: 


SS 2 


d ðu 
4 6,u2)=2 | (22) a 
aa 3U”) ay Mg 
0 
care reprezinta ecuatia integrala a energiei pentru stratul limita bi-dimensional, laminar in curgere 


incompresibila. 
In caz de curgere turbulenta, ecuatia integrala a energiei va avea forma 
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0O 
d To Ou 
— (63U3) = 2 | —— d 
0 
Pentru a vizualiza grosimea de deplasare, grosimea de impuls si grosimea de disipare a energiei, este 
convenient sa se reprezinte pentru cazuri simple de distributie liniara a vitezei, ca in figura de mai jos. 


In figura este prezentat un strat limita cu distributie de viteza liniara in care: 
ô — grosimea stratului limita 


1 

6, = 5 6 este grosimea de deplasare 
1 l 

62 = = 6 este grosimea de impuls 


1 l l 
Ô, = —6 este grosimea de energie 
4 
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Desprinderea stratului limita la curgerea in jurul unui corp 


Experimental se stabileste usor ca in spatele unui corp exista o zona cu un pronuntat caracter 
rotational, ceea ce modifica distributia vitezelor in fluidul liber fata de cazul cand aceasta nu ar exista. 
Se spune ca are loc desprinderea sau separarea stratului limita de suprafata corpului. 


Fie un corp de sectiune eliptica in curent uniform la mari distante amonte, paralele cu axa AA’ a 
elipsei. Datorita simetriei vom considera doar miscarea din semiplanul superior lui AA’. Daca fluidul 
ar fi perfect, viteza tangentiala la conturul ABA’ mai intai va creste de la zero in punctul A pana la 
valoarea maxima in punctul B, cel mai indepartat de axa mare, ca apoi sa scada pana la zero in A’. 
Presiunea va scadea deci cand particula fluida se deplaseaza din A spre B si apoi va creste cand 
particula se deplaseaza mai departe din B in A’. 

e Punctul A — de extremitate a bordului de atac unde presiunea statica este maxima — este 
punctul in care practic aproape toata energia cinetica de deplasare a fluidului din linia de curent 
perpendiculara pe corp este transformata in energie statica de presiune 

e Portiunea AB — zona de curgere laminara vascoasa 

Daca fluidul este vascos atunci distributia presiunii in stratul limita e practic egala cu distributia 
presiunii in miscarea fluidului liber. In ce priveste distributia vitezei, ea va fi esential modificata. O 
particula care se gaseste in stratul limita va fi franata in miscare datorita fortelor de frecare. 

e Portiunea BC — zona de scadere a adeziunii stratului limita la perete 

Incepand din punctul B coeziunea in interiorul stratului limita scade, insa acesta ramane stabil pana 
in punctul C unde incepe ruperea acestuia. 

e Punctul C — de desprindere (rupere) a stratului limita - gradientul advers de presiune impus 
de potentialul de curgere exterior determina ruperea stratului limita si aparitia unei curgeri in 
sens invers. 

e Portiunea CA’ — zona de desprindere si indepartare a stratului limita de perete 

Cum in portiunea BA’ atat fortele de frecare cat si presiunea se opun miscarii, urmeaza ca viteza se va 
anula inainte ca particula sa ajunga in A’, deci fie acel punct C. Intre C si A’ sensul miscarii se 
schimba, astfel incat linia pe care avem viteza nula, ca si pe suprafata corpului, este o curba CC’ care 
se intinde in avalul corpului. In punctul C deci are loc desprinderea stratului limita. 
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De obicei punctul de desprindere se afla in punctul pana la care e posibil calculul stratului limita. Din 


cele aratate, acest punct va fi acela in care 
( ) = 0 
dy 


y=0 
deci in care efortul tangential t = 0. 


In figura de mai jos este aratat comportamentul liniilor de curent in apropierea punctului de 
desprindere, precum si variatiile grosimii stratului limita in apropierea acestui punct. 


i 


eh EA 


\ 


| 
AN 
NANNY 
IN 
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HI 
| 
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| 
N 


E de remarcat cresterea notabila a grosimii stratului limita in apropierea punctului C, astfel incat 
mai departe ecuatiile folosite pana acolo, bazate pe ipoteza ca 6 « 1 isi pierd valabilitatea. 

Un fapt demonstrat, cum ca desprinderea are loc in regiunea miscarii intarziate dp/dx>0, se poate 
deduce si pe alta cale. 

Pentru miscarea stationara, deoarece u(x,0)=v(x,0)=0, urmeaza o forma echivalenta: 


(=) dp 
v| — = — 
dy? a dx 


Deci in apropierea suprafetei corpului, curbura profilului vitezelor este determinata numai de 


presiune. 
Analizand si variatia raportului dp/dx apar urmartoarele situatii: 


«0? . n 
e In regiunea unde dp/dx<0 avem si E < 0 peste tot in stratul limita. 


; . (8? d 
e In regiunea unde dp/dx>0, deci (=) > 0, cum la o anumita distanta de perete intotdeauna 
y=0 

07u i . l 07u . i . i 
3y < 0 va trebui sa existe un punct in care a O ; profilul vitezelor va avea in acest caz in 
stratul limita un punct de inflexiune. In punctul de desprindere tangenta la profilul vitezelor 
face un unghi nul cu normala la perete, deci trebuie sa fie neaparat punctul de inflexiune al 
acestui profil. Ecuatia impulsului unei linii de curent din interiorul stratului limita este 


aproximativ dupa ecuatia: 


Unde s, y sunt coordonatele liniei de curent respectiv coordonata normala la suprafata. Astfel, acest 
d ae 
gradient advers de presiune este cand Pa > 0, care poate cauza ca viteza u sa scada de-a lungul lui s si 


posibil sa ajunga la zero daca gradientul de presiune este suficient de puternic. 


ðu 3? u 
dy’ ðy2 ` 


In urmatoarea figura de mai jos sunt reprezentate variatiile lui u, 


Pentru corpul din figura experienta arata ca desprinderea vartejurilor are loc alternativ, odata in C, 
apoi in D, pe urma in C, s.a.md., astfel incat in avalul corpului apare o configuratie a vartejurilor in 
care vartejurile care s-au desprins din C se gasesc la inceput la distante practic egale de vartejurile 
vecine, ca si vartejurile care s-au desprins din D, iar dreapta ce trece prin centrul unui vartej desprins 
din D si este perpendiculara pe AA’ se gaseste intre doua vartejuri desprinse din C, la inceput practic 
la distanta egala de acestea. Configuratia mentionata (“in sah”) se numeste “aleea lui Karman”, si 
studii asupra e1 in hidrodinamica fluidelor perfecte au fost intreprinse de multi cercetatori, printre care 
mentionam pe Blagovest Dolapciev. 

Fenomenul de desprindere al vartejurilor cel putin pentru miscarea plana, este in mod esential 
nestationar, ceea ce complica si mai mult abordarea anumitor probleme. 

Pentru indepartarea desprinderii sau pentru deplasarea punctului C inspre avalul curentului s-au 
propus diferite metode, printre care si miscarea pe o placa plana. 
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Exemplu de simulare numerica a curgerii bidimensionale in jurul unei 
elipse 


Date de intrare: 
e Corp eliptic: Lxl=2m x 1m 
e Viteza in sectiunea de intrare: v=l m/s 


Meshatura sistemului studiat 


Fiind vorba de o geometrie neliniara a corpului studiat in cazul de fata este recomandat un mesh cu elemente 
triunghiulare. 

Spatiul de meshatura ocupat de fluid in curgere in jurul corpului este luat suficient de mare, incat marginile 
exterioare ce il delimiteaza sa nu influenteze semnificativ modelul. Bineinteles, densitatea de elemente finite va 
creste odata cu apropierea de peretele corpului, pentru a se asigura o retea cat mai fina de noduri in zona 
stratului limita. 
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Campul de viteze al fluidului din jurul corpului elipsoidal 


Observand rezultatele obtinute pentru distributia campului de viteze, se disting 3 zone importante: 


e Zona bordului de atac — caracterizata de o presiune statica crescuta rezultata din conversia energiei 
cinetice a fluidului ce intalneste corpul in energie potentiala din forte de presiune. Aceasta presiune 
statica rezulta din procesul de franare a fluidului in zona vecinatatii corpului. 

e Zona laterala intermediara a corpului — specifica formarii si evolutiei stratului limita 

e Zona bordului de fuga — caracterizata de turbulente si este rezultata din indepartarea stratului limita de 
corp. 
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Campul de presiuni in jurul corpului elipsoidal 


Este de remarcat cum presiunea atinge un maxim in zona bordului de atac, ulterior urmand o descrestere 
acesteia pana in punctul critic din care are loc ruperea stratului limita. 
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